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APPARATUS  FOR  MEASUREMENTS  OF  TIME  AND  SPACE  CORRELATION*! 
By  A.  Favre,  J.  Gavlglio,  and  R.  Dumas 

SUMMARY 


These  researches  are  made  at  the  Laboratoire  de  Mecanique  de 
1 'Atmosphere  de  I'l.M.F.M.,  for  the  Office  National  d'Etudes  et  de 
Recherches  Aeronautiques  (O.N.E.R.A.) ^  with  the  aid  of  the  Ministere  de 
I'Air,  and  of  the  Centre  National  de  la  Recherche  Scientifique  (C.N.R.S.)» 

We  have  alxeady  explained  in  previous  communications  (refs.  1  to  22) 
the  principle  and  have  given  description  of  the  apparatus  for  statisti- 
cal measurements  of  time  and  space  correlation  between  two  random  vari- 
ables considered  with  a  relative  delay  T  in  time,  which  we  have  designed 
with  a  special  view  to  the  study  of  tiorbulence. 

We  shall  briefly  review  now  the  improvements  brought  to  the  experi- 
mental apparatus,  the  control  of  the  measurements,  and  a  few  particulaj* 
applications . 


EXPERIMENTAL  APPLIANCES 
Wind  Tunnel  (Refs.  2,  5,  ^4-,  6,  12,  and  15) 


The  low- turbulence  wind  tunnel.  Installed  indoors,  has  been  trans- 
formed (fig.  1).  The  test  section  A  is  80  x  80  X  270  cm;  the  cone  B 
with  low  gradients  has  a  contraction  ratio  of  15-  It  is  preceded  by  a 
settling  chamber  C,  in  which  are  set  five  removable  damping  screens, 
and  by  a  rapid  expansion  diffuser  D  with  a  screen  analogous  to  those  of 
the  low-tvirbulence  wind  tunnel  of  the  R.A.E.  (ref.  25).  The  intensity 
of  turbulence  is  of  the  order  of  0.00210  with  the  diffuser  screen  and 

*"Appareil  de  Mesures  de  la  Correlation  Dans  le  Temps  et  L'Espace." 
La  Recherche  Aeronautique  No.  51,  Jan. -Feb.  1955,  PP-  57-^' 

-'-A  Commiinication  to  the  8th  International  Congress  on  Theoretical 
and  Applied  Mechanics,  Istanbul,  Aug.  1952. 
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of  0.00150,  0.00090,  0.00055,  0.00014-5,  and  O.OOO58  with  1,  2,  5,  ^, 
and  5  damping  screens,  respectively. 

At  the  entrance  of  the  test  section  may  be  placed  one  of  the  two 
tiorbulence  grids  similar  to  those  of  the  NBS  with  elements  of  O.626  inch 
and  0.196  inch  in  diameter  and  with  meshes  M  of  5  lA  inches  and  1  inch, 
respectively. 

A  flat  plate  of  1,9^0  x  80  cm  can  be  set  in  the  test  section  for 
the  study  of  the  boiindary  layer . 

Hot-Wire  Anemometers  (Ref s .  2,  3,  5,  7,  8,  11,  I7,  and  20) 

Two  hot  wires  can  be  moved  along  the  axis  X^     of  the  test  section 
and  along  the  orthogonal  axis  X5 . 

The  amplifiers  reproduce  the  oscillations  with  an  approximation  of 
±6  percent  in  energy  and  for  the  band  pass  1  to  6,000  cps,  more  recently, 
for  1  to  3,500  cps  with  tk   percent.  The  compensation  is  made  by  the 
sqiiare  waves  method. 

Figure  2  shows  the  two  amplifiers  A,  the  power  supplies  B,  the 
square  waves  generator  C,  and  a  generator  of  sine  waves  from  0.6  to 
30,000  cps  D. 


Apparatus  for  Statistical  Measurements  of  Time  and 

Space  Correlation  (Refs.  1,  2,  3,  5,  8,  9,  10,  11,  and  I3) 

The  apparatus  records  simultaneously  and  reproduces  with  a  relative 
delay  T  the  two  voltages  coming,  for  instance,  from  the  two  anemometers 
and  measures  the  time- correlation  coefficient  (refs.  1,  2,  and  5). 

It  must  be  borne  in  mind  that,  for  this  measurement,  it  is  necessary 
and  sufficient  that  the  phase  shifts  of  the  component  waves  due  to  the 
inevitable  transformations  should  be  identical  in  the  two  channels  for 
each  frequency;  the  difference  T  of  the  delays  systematically  intro- 
duced being  the  same  for  all  the  frequencies . 

The  voltages  are  recorded  on  magnetic  tape: 

(a)  directly  for  the  waves  of  frequencies  comprised  between  200 
and  2,500  cps 
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(b)  on  an  amplitude  modulated  wave  of  5,500  cps  frequency,  demod\ila- 
ted  after  reproduction,  for  the  waves  of  frequencies  comprised 
between  1  to  25  cps 

(c)  by  overlapping  of  the  two  devices  from  25  to  200  cps. 

The  band  pass  is  therefore  unbroken  from  1  to  2,500  cps  (refs.  8, 

9,  10,  and  11)  , 

The  apparatus  consists  of  (fig.  2): 

The  two  recording  preamplifiers  E,  with  a  common  oscillB.tor  of 
3,500  cps,  and  on  either  channel  the  networks  effecting  the  division  of 
the  respective  spectra  into  waves  of  200  to  2,500  cps  for  direct  recording, 
of  1  to  25  cps  for  modulation  of  the  5>500  cps  carrier,  and  of  25  to 
200  cps  for  the  aforementioned  overlapping 

The  two  recording  amplifiers  F  (with  appended  oscillators  of 
80,000  cps  and  of  i+0,000  cps)  and  their  power  supplies  G 

The  two  magnetic  recorders-reproducers  H,  one  of  which  has  a 
device  I  for  adjusting  the  length  of  the  tape  between  the  recording 
and  play-back  heads  and,  consequently,  the  delay  T  (refs.  2  and  5) 

The  two  play-back  preamplifiers  F 

The  two  play-back  amplifiers  J,  with  their  networks  effecting, 
after  the  detection  of  the  carrier,  the  reconstitution  of  the  initial 
spectrum  on  each  channel 

The  two  power  amplifiers  K  which,  after  adding  on  one  hand  and 
subtracting  on  the  other,  supply  the  thermocouples;  the  output  voltages 
of  these  thermocouples  represent  the  mean  squares  M^^  of  the  sum  and 
Mh  of  the  difference,  respectively  (Dryden  method);  two  extra  amplifiers 
for  the  oscilloscope 

The  apparatus  L,  which  gives  the  sum  and  the  difference  of  M^ 
and  M-)-,  (and  will  soon  reckon  their  quotient),  together  with  an  oscil- 
lator for  20,  200,  1,000  cps,  by  means  of  which  may  be  effected  a  sup- 
plementary control  of  the  phase  and  energy  response  of  the  two  channels 
during  the  process  of  measurements;  this  oscillator  is  also  used  for 
the  measurement  of  the  delay  T 

The  power  supplies  N;  the  oscillograph  0;  the  galvanometer  M 
with  adjustable  period  (5  sec). 

The  equality  of  phases  on  the  two  channels  is  checked  for  each 
frequency  (refs.  2  and  3) • 
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When  the  two  voltages  correspond  to  two  variables  measured  at  two 
different  points  of  space,  one  obtains  the  coefficient  R  of  time  and 
space  correlation. 

In  the  particular  case  when  it  represents  the  same  variable  measured 
at  the  same  point  of  space,  one  gets  the  autocorrelation  coefficient 
from  which  the  spectral  function  can  be  deduced. 

Spectral  Analyser  (Ref s .  9,  10,  15,  Ik,   15,  and  l6) 

A  selective  feed-back  ajnplifier  has  been  designed  in  order  to  deter- 
mine the  spectral  function  directly,  apart  from  the  autocorrelation 
method. 

The  band-pass  is  in  width  4.5  percent  of  the  median  frequency 
chosen  for  frequencies  comprised  between  I.5  and  2,500  cps. 

The  energy  is  measured  by  means  of  a  thermocouple  and  a  long-period 
galvanometer  (19  sec). 

COIWROL  OF  THE  MEASUREMENTS  (REFS.  2  to  I7,  AND  21) 


I.  The  response  curve  of  the  amplifiers  of  the  hot-wire  anemometers 
(ref.  7)  has  been  determined  as  regards  energy  and  phase  by  means  of 
sine  waves  and  square  waves  of  variable  frequencies,  and  so  has  the 
band-pass  of  the  spectral  analyser. 

II.  The  total  response  curve  of  the  apparatus  for  measurements  of 
time  correlation,  together  with  the  hot-wire  anemometers,  has  been  estab- 
lished with  sine  waves  of  variable  frequencies,  as  regards  the  energy 
(fig.  5),  the  equality  of  phases  on  the  two  circuits  for  each  elementary 
oscillation  (refs.  2  and  5)  (the  delay  T  being  a  multiple  of  the  period), 
and  the  autocorrelation  function  (refs.  2  to  15)- 

The  energy  is  constant,  variations  not  exceeding  ±l6  or  19  percent 
of  the  standard  deviation  0.05  for  1  to  2,500  cps  (refs.  2  to  16) . 

Frequent  checkings  are  made  in  the  course  of  measurements  especially 
in  the  case  of  20,  200,  and  1,000  cps  oscillations  frequencies. 

III.  The  spectra  are  obtained  on  one  hand  by  Fourier  transform  of 
the  autocorrelation  curves  and  on  the  other  hand  directly  by  means  of 
the  spectral  analyser  and  are  then  systematically  compared.  These 
measurements  are  made  in  the  case  of  turbulence  downstream  of  a  grid 
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(refs.  9,  10,  15,  1^,  and  2h)   and  also  in  the  case  of  a  boundary  layer 
of  a  flat  plate  (examples:   figs.  ^4-,  5,  and  8)  (refs.  I5  and  2k) 


e  spectra  are  satisfactorily  concordant  from  5  "to  2,000  cps 
4). 


Thes 

now  (fig.  4) 

rV.  The  vibrations  of  the  hot-wire  holders  (ref .  2)  and  of  the 
tiinnel  walls  (ref.  W)   have  been  examined  and  reduced;  their  effects  are 
watched  with  the  oscilloscope  dioring  the  experiments. 

The  effect  on  the  measiorements  of  the  acoustic  resonances  has  been 
found  to  be  negligible  (refs.  6  and  12). 

V.  Comparative  measurements  have  been  made  (refs.  \,    6,  9,  10,  12, 
and  15)  at  each  stage  of  the  transformations  of  the  wind  tunnel: 

(a)  Length  of  the  test  section  reduced  from  400  to  270  cm, 
contraction  ratio  of  the  cone  increased  from  GX   to  I5, 
the  new  settling  chamber  placed  behind  the  diffuser  with 
its  screen,  and  five  damping  screens  inserted  successively 
in  the  bulge 

(b)  Installation  of  the  first  paj-t  of  the  retxirn  circuit  with 
two  corners  with  vanes 

(c)  Size  of  the  test  section  reduced  from  80  x  80  cm  to 
80  X  40  cm  (ref.  I5) . 

These  comparisons  have  not  shown  any  appreciable  changes  in  the 
results  from  one  stage  to  the  other,  as  regards  the  turbulence  behind 
a  grid. 

One  must  except,  however,  the  influence  of  the  initial  turbulence 
and,  consequently,  of  the  damping  screens,  which  is  very  important  in 
the  case  of  the  boundary  layer  and  the  shape  of  the  grid. 

VI.  The  comparison  of  the  time  correlation  curves,  obtained  with 
different  hot  wires,  disclosed  appreciable  differences  for  wires  of  the 
same  diameter  and  especially  for  wires  of  different  diajneters .  The 
improvement  of  the  compensation  by  the  squaire-waves  method  has  made  it 
possible  to  reduce  those  differences  (fig.  6)  (refs.  10  and  13). 

Aerodynamic  taring  attempts  of  the  hot-wire  anemomenter  are  in 
progress,  the  wire  being  placed  in  the  field  of  a  propeller  (refs.  I7 
and  20) . 

The  measurement  errors  seem  to  be  mainly  due  to  the  hot  wires;  the 
comparative  experiments  are  preferably  made  with  practically  identical 
wires . 
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VII.  The  results  obtained  for  the  turbulence  behind  a  grid  with  the 
apparatus  for  time  and  space  correlation  measurements : 

transversely,  i.e.  at  several  distances  Xj  (refs.  9>  10>  I6, 
21,  and  2k) 

longitudinally,  i.e.  at  several  distances  Xi  (refs.  1^,  I6, 
and  21) 

the  delay  T  being  zero  and  agreeing  with  the  results  con- 
cerning the  transverse  and  longitudinal  space  correlations  ^ 
and  f,  respectively,  (example:  fig.  7)  (refs.  9^  10^  21, 
and  2k). 

VIII.  The  comparison  of  the  spectra  (fig.  8)  (refs.  I5  and  2k)   with 
those  of  H.  L.  Dryden  (NBS)  and  of  F.  G.  Simmons  (NPL)  established  for 
the  t\irbulence  behind  a  grid  in  quite  similar  conditions  is  satisfactory. 
The  measurements  have  been  especially  developed  up  to  I.5  cps .  The 
comparison  with  A.  A.  Townsend's  spectra  (ref .  22)  is  also  satisfactory 
although  the  experimental  conditions  are  slightly  different. 

IX.  The  qualitative  comparison  of  the  spectra  established  for  the 
turbulence  of  a  flat-plate  boiindary  layer  (fig.  9)  (ref.  I5)  with  those 
of  A.  A.  Townsend  is  fairly  satisfactory,  although  the  experimental 
conditions  present  some  differences  (ref.  25). 


A  FEW  APPLICATIONS  OF  TIME  CORRELATION  MEASUREMENTS 


Among  the  applications  which  can  be  made  of  time- correlation  measiore- 
ments  between  two  variables  or  of  autocorrelation  for  one  variable,  and, 
in  addition,  to  those  concerning  the  turbulence  studied  elsewhere  (ref.  22), 
we  shall  mention  some  other  experiments  (refs.  I8  and  19)- 


Application  of  Time  Correlation  to  the  Relative  Increase  of  a 

Periodic  Signal  Disturbed  by  a  Random  Parasite  (refs.  3,  18,  and  19) 

1.  In  a  previous  paper  (ref.  3)  we  have  considered  the  coefficient 
of  time  correlation  R]^]^  22  between  two  signals  uj^  and  U2  disturbed 
by  two  parasites  e^  and  62,  respectively,  and  the  coefficient  of  time 
correlation  Rj^q  20  between  the  nondisturbed  signals^  and  noticed  that, 
if  one  can  choose  a  delay  T,  such  as  the  time  correlation  coefficient 
between  the  parasites  and  between  the  signals  and  the  parasites,  to  be 
negligible 


With  ui  =  U2  =  e2_  =  £2  =  0,   q-^  =  ^±/^±^}      '^2  -   ^2^ 


^2 
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u^eg  =  e 


1U2  =  u-^e-^   =  U2e2  =  £162  =  0  (l) 


one  finds  that 


Rll,22  =  1^10,20  lAl^/l/(^  +  I1)  (^  +  ^12) 


The  coefficient  of  time  correlation  between  the  disturbed  signals 
is  proportional  to  the  coefficient  of  time  correlation  between  the  non- 
disturbed  signals,  when  the  delay  is  such  that  the  coefficients  of  time 
correlation  between  the  parasites  and  between  the  signals  and  the  parasites 
are  zero . 

With  periodic  signals  and  random  parasites,  condition  (l)  is  easily 
fulfilled;  thus,  these  signals  can  be  detected.  With  sine  waves  of  the 
same  frequency,  the  coefficient  of  time  correlation  R]_i  22  °r  o^  auto- 
correlation Rxi  11     will  be  for  appropriate  delays  a  sine  wave  of  the 

same  frequency,  which  will  constitute  a  time  detection  of  these  signals 
(ref.  18). 

In  the  case  of  a  single  signal,  it  may  be  preferable  to  replace  the 
autocorrelation  Rix  ]^]_  by  the  time  correlation  B.-^^^   12  with  an  iden- 
tical signal  disturbed  by  a  less  intense  parasite. 

2.  We  have  made  trials  of  the  detection  of  a  sine  wave  by  autocor- 
relation (refs.  18  and  19).  For  example,  figure  10  (ref.  I9)  shows  the 
autocorrelation  curves  of  the  parasite  alone  —  the  spectrum  of  which 
stretches  mainly  from  UOO  to  1,100  cps  —  and  of  the  signal  of  ~  6OO  cps 
distiirbed  by  the  parasite . 

It  may  be  seen  that,  for  a  delay  T  superior  to  8  ms,  the  latter 
is  close  to  the  sine  wave  of  ~  60O  cps  of  amplitude  O.k. 

Other  tests  have  enabled  us  to  verify  the  above  formula  and  to 
detect  the  presence  of  the  signal  for  very  small  values  of  q  with  a 
great  sensitivity  (ref.  19) . 
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Reduction  of  Parasites  Ratio  With  Reference  to 

Disturbed  Phenomena  by  Addition  With  Delay  (Ref s .  18  and  I9) 

1.  Let  ut  be  the  variable  representing  a  phenomenon  disturbed  by 
a  parasite  e-^.     Let  us  add  the  disturbed  variable  to  itself  with  a 
time  delay  T;  we  then  have^ 

Q  ^  ^  +  %0,10 

^     1  +  f'oi^oi 

Thus,  it  is  possible  to  reduce  the  parasite  disturbing  a  phenomenon 
by  addition  with  delay.  This  reduction  will  be  at  its  maxim\im  when  T 
is  such  that  the  ratio  of  the  autocorrelation  coefficients  of  the  phenom- 
enon, and  of  the  parasite,  to  both  of  which  the  unity  has  been  added, 
will  be  at  its  maximum. 

In  the  case  of  a  periodic  variable  u-t  and  of  a  random  parasite, 
one  may  take  a  delay  multiple  of  the  period  and  such  as 

Roi,oi«o 

then, 

Q, 

-  »2 

<1 

It  may  be  preferable  to  filter  the  parasite  so  that,  for  a  T 
multiple  of  the  period,  one  should  have 


^01,01 
then 


'■> 


1 

The  operation  may  be  repeated  in  series,  choosing  the  optimum 
delays;  one  may  also  make  several  additions  in  parallel  with  optimum 
delays . 

\ith  £^2  =  et+T^,   ut2  =  ut+T2,   q  =  u^^U^^ , 


Q  =  (ut  +  n^+^f/[e^  +   et+T)2 
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2.  We  have  made  (ref.  18)  some  measurements,  using  a  sine  signal  of 
frequency  n  and,  for  the  parasite,  the  above-mentioned  tiorbulent  veloc- 
ity fluctuations . 


For  a  delay  T  such  that 


^01,01  -  0 


and 


^10,10  ~  1 

n 157  180  210  500  570  570  570 

Q/q 1.90  1.97  2.0^  2.41  2.0i+  I.9I  2.22 

the  mean  value  2.08  approximates  the  theoretical  value  2. 

For  a  delay  T  such  that 

1^01,01  ~  -0-5 

and 

^10,10  ~  1 

and  n  =  80,  the  ratio  Q/q  reached  4.25,  4.52,  and  4.42,  slightly 
above  the  theoretical  value  4. 

Translated  by  A-  Favre 
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Figure  1.-    Low -turbulence  wind  tunnel. 


Figure  2,-    Apparatus  for  measurements  of  time  correlation. 
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Figure  6.-    Autocorrelation  measured  by  different  hot  wires  and  by 
two  methods  of  compensation. 
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Figure  7.-    Control  of  time  and  space  correlation  measurements  with 
zero  delay,  transversely  g,    longitudinally   f ,  downstream  of  a  grid. 
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(c)     IMA  square  vave  compensation,      (d)     j^  semlempirical  ccoipensation. 


Figure  8,-    Comparison  of  turbulence  spectra  downstream  of  a  grid 
obtained  at  NBS,  NPL,  and  LMA. 
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